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ABSTRACT 

The progenitor stars of several Type Ilb supernovae (SNe) show indications for extended hydrogen envelopes. 
These envelopes might be the outcome of luminous energetic pre-explosion events, so-called precursor erup¬ 
tions. We use the Palomar Transient Factory (PTE) pre-explosion observations of a sample of 27 nearby Type 
lib SNe to look for such precursors during the final years prior to the SN explosion. No precursors are found 
when combining the observations in 15-day bins, and we calculate the absolute-magnitude-dependent upper 
limit on the precursor rate. At the 90% confidence level. Type lib SNe have on average < 0.86 precursors as 
bright as absolute R-band magnitude —14 in the final 3.5 years before the explosion and < 0.56 events over the 
final year. In contrast, precursors among SNe Iln have a > 5 times higher rate. The kinetic energy required to 
unbind a low-mass stellar envelope is comparable to the radiated energy of a few-weeks-long precursor which 
would be detectable for the closest SNe in our sample. Therefore, mass ejections, if they are common in such 
SNe, are radiatively inefficient or have durations longer than months. Indeed, when using 60-day bins a faint 
precursor candidate is detected prior to SN 2012cs 2% false-alarm probability). We also report the detection 

of the progenitor of SN 201 Idh which does not show detectable variability over the final two years before the 
explosion. The suggested progenitor of SN 2012P is still present, and hence is likely a compact star cluster, or 
an unrelated object. 

Subject headings: stars: mass-loss — supernovae: general — supernovae: individual: SN201 Idh, SN2012P, 
SN2013bb, SN2012CS 


1. INTRODUCTION 

Violent stellar activity, as for example precursor eruptions, 
during the years prior to the terminal explosion of a star 
as a supernova (SN) can give insights into the properties of 
progenitor stars at this crucial stage of their life. In case 
these eruptions are causally connected to the final explosion 
they might even provide hints about the explosion mecha- 


viously, mainly among Type Iln SNe (e.g. , 
Pastordlo et al.l l2007t [Eraser et alJ 120131: 


Folev et al.l 

2007 

Ofek et al.l 

2013 


also reported for a Type Ic SN (ICorsi et al.ll2014li . Here 


was 

we 
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search for precursor eruptions, i.e., small explosions prior to 
the SN, in the pre-explosion light curves of Type lib SN. 

The defining property o f SNe lib, the firs t known example 
of which was SN 1987K (lFilippenkolll988h . is their spectral 
evolution. The spectrum first resembles that of a hydrogen- 
rich Type II SN and later turns into a helium -dominated spec¬ 
trum simila r to that of a Type Ib SN (e.g., iFilippe nko et alJ 
1993L11994: IChornock et al.ll201 lb l^n-Ami et alJl2015t see 


Filippenl^ 19971 for a general review of SN spectra). This 


evidence suggests that the progenitors of SNe lib likely con¬ 
sist of a helium core surrounded by a low-mass hydrogen 


tures fade awav within a few weeks (e.g.. iWooslev et al. 

19941 iTsvetkov et al.l I2009t IChevalier & Soderberg 

201C; 

Bersten et al. 12012 

lErgon et alJl2014allbt iMarion et all 

2014; 

Nakar & Piroll2014 

). 


Another characteristic is that some SNe lib show 
two optica l photometric peaks in their light curve (e.g., 
SN 19 93J, I Wheeler et all [19^ SN 2 01 Idh, lArcavi et alJ 
IWIt SN 2013df. IVan Dvk et al.llMl) . The first peak can 
be attributed to thermal e mission from a shock-heated ex¬ 
tende d hydrogen envelope (iBersten et al.ll201^ iNakar & Pirol 
l2014h . whil e the second one is presumably powered by ra¬ 
dioactivity. iNakar & Pirol (1201 4l) suggest that the envelope’s 
radius and mass can be inferred from the shape of the first 
peak. They have also shown that standard red giant or Wolf- 
Rayet progenitors cannot explain the observed light curves. 
Indeed, large photospheric radii of several hundred solar 
radii have been measured for several SNe Hb at early times 

(e.g., |^)oslev et al.||1994lTsvetkov et al.||20( )9t IBersten et alJ 

l2012l:[Horesh et alJl2013HErgon et al.ll2014allbli. 

Additionally, an ultraviolet (UV) excess relative to the 
spectral shape in the optical has been observed for several 
SNe lib whose UV spectra were obtained within the first 
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month after explosion. It can be explained by the presence 
of a de nse circumstellar med ium (CSM) above the photo¬ 
sphere (iBen-Ami et al.llT015l) . Double-peaked light curves 
and exces s UV radiation are , however, not observ ed for all 
SNe Ilb (lArca vi et alJ 12015: iB en-Ami et al.l 12015]) . More¬ 
over, IChevalier & Soderberg (l2010l) argue that some SNe lib 
show indications for much more compact progenitors having 
R^IRq. 

It is unclear how extended envelopes or dense CSM shells 
are created. They cou ld consist either of mat e rial in hydro¬ 
dynam ic equilibrium (iBenvenuto et al.l 12013l iBersten et alJ 
12012h or of unbound material ejected during th e interaction 
with a binary companion (e.g.. IChevalieiil2012h . by a stellar 
wind, or by an in stability related to the final stages of stel¬ 
lar bur ni ng (e.g., [Rakavv & Shavivlll967t lArnett & MeakinI 
l2011allbl : IOuataert & Shiodell2012HShiode & Ouataerlll2014ft . 

Such an instability can generate a precursor event that can re¬ 
lease some of its energy in visible light. 

A systematic search for precurs or eruptions amon g SNe Iln 
has previously been performed bv lOfek et al.l (l2014t) . with six 
likely precursor explosions observed in a sample of 16 ob¬ 
jects. A quantitative analysis showed that, assuming a uni¬ 
form population, at the 99% confidence level at least 98% of 
the SNe Iln exhibit one or several precursor events brighter 
than an absolute magnitude of —14 within the last 2.5 years 
prior to their explosion. These results apply to all SNe Iln 
under the assumption that they form a homogeneous popula¬ 
tion w hich is wel l repre sented by the used sample. Further¬ 
more, lOfek et akl (120141) find possible correlations between 
the amount of energy radiated in a precursor event and the SN 
rise time, peak magnitude, and radiated energy. If real, the 
correlation is consistent with the idea that SNe Iln are pow¬ 
ered mainly by the interaction of the SN ejecta with a massive 
CSM created in precursor eruptions. 

Only SNe having strong and pe rsistent narrow hy drogen 
lines were included in the sample of lOfek et al.l (I2Q14I) . It was 
found recently that some core-collapse SNe, dubbed ’’flash- 
spectroscopy SNe”, exhibit narrow hydrogen lines within 
the first days after the explosion which va nish subsequently 
(iGal-Yam et al.l l2Q14t iKhazov et al.l l2015h . This evolution 
might indicate that a dense CSM shell or an envelope is lo¬ 
cated only in the immediate vicinity of the progenitor star and 
is swept up by the ejecta shortly after the explosion. Owing 
to the short duration of the CSM interactio n, these objects are 
not considered in the sample of SNe Iln of lOfek et al.l (l2014t) 
and might not undergo precursor explosions as frequently. We 
note that the detection of flash-spectroscopy signatures pro¬ 
vide further evidence for mass ejections prior to SN explo¬ 
sions. 

In another search for precursors among SNe Iln by 
iBilinski et al.l (12015h . no pre-explosion events were found. 
However, with six objects their sample is smaller, and they 
do not provide an absolute-magnitude-dependent upper limit 
on the precursor rates, so a close comparison is not possible 
at this stage. 

Here, we extend the search of lOfek et al.l ([203) to Type 
lib SNe. In Sj2]we explain how the SN sample was selected, 
and in 0we describe the observations used in this analysis. 
Our search for precursor explosions is presented in 0 
The pre-explosion observations of three nearby SNe are 
evaluated in more detail in 0 In 0 we calculate the sample 
control time and derive an upper limit on the precursor rate 
of SNe Ilb. Section |7] discusses whether the ejection of a 
low-mass stellar envelope is likely detectable in this search. 


and ij^summarizes the results. 


2. SAMPLE SELECTION 


The sample is selected by initially considering all nearby 
SNe I lb detected either by the Palom ar Transient Eactory 
(PTE: iLaw et al.ll2009l: iRau et al.ll2009^ or announced in As¬ 
tronomer’s Telegrams or lAU circulars since 2009 (i.e., after 
the start of the PTE project). SNe Ilb are chosen based on 
their spectral evolution, and the main criterion is the appear¬ 


ance of helium features at 5876 A, 6678 A, and 7065 A a 
few weeks after the explosi on. Tools such as the S upernova 
Identi fication Code (SNIP ; [Blondin & Tonrvll2007l) and Su- 
perfit (iHowell et al.ll200-^ are used to compare our spectra to 
the spectra of known SNe Ilb in case the spectral features are 
not obvious. We note that in general, our SN classification is 
based on human decision and hence may be biased. To in¬ 
crease the chances of detecting faint precursors, we restrict 
ourselves to nearby SNe with z ^ 0.05, corresponding to a 
luminosity distance of 220 Mpc. Only SNe having a large 
number of PTE observations prior to their explosion are se¬ 
lected. We require about 20 images (either before or long 
after the SN explosion) to construct a high-quality reference 
image, and another 20 science images before the explosion 
date. 

Our final sample consists of 27 SNe listed in Table [T] The 
majority were discovered by PTE, and additional SNe were 
found by the Lick Observatory Sup ernova Search (LOSS; 
iLi et alJ [20001: lEilinnenko et al.ll200ih . the Puckett Observa- 
tory Supernova Se arclF^ the La Silla Quest survey (LSQ; 


iBaltav et~al 


l2013h. the Catali na Real-Time Transient Sur¬ 


vey (CRTS; Drake et ani2009al) . the Italian Supernova Search 
ProiecPl (ISSP), and the Panoramic S urvey Telescope an d 
Rapid Response System (Pan-STARRS: fHodapp et alj|2004l) . 
Eor three SNe (SN2011ef, SN2012an, and LSQ 12fwb), no 
spectra are publicly available and we rely on the classifica¬ 
tion published in Astrono mer’s Telegrams and Ce n tral Bureau 
Electr onic Telegrams bvlB anchard et'H] (1201 ll) . IChen et ^ 
(|20I3), and lHadiivska et ^ ( 20121) . Representative spectra of 
all SNe discovered by PTE are shown in Eigures IMTOl and 


are electronically availab le from the WISeREP webpagj^ 
(lYaron & Gal-Yarnll2012^ . The spectra were acquired using 
various facilities which are listed in Tabled 
Our sample includes SN2013cu (PTF13a st), which has 
been classified as a flash -spectroscopy SN (iGal-Yam et alJ 
l2014tl^azov et aTll20L5l) . After the initial CSM interaction, 
SN2013cu evolves into a SN Ilb (see the spectrum in 
Figure fTOl l and is thus included in this sample. Out of all 
core-collapse SNe having spectra within the first ten days 
after their l ast nondetection, ^ 14% show flash-spectroscopy 
signatures (IKhazov et al.l I20f5h . It is therefore likely that 
other SNe in our sample for which no early spectra have been 
obtained belong to this group. 


3. OBSERVATIONS 

The observations used here were obtained with the 48-inch 
Oschin Schmidt telescope at Palomar Observatory (P48), as 
part of the PTE project. PTE searches for transient sources 

Fittp://www.cometwatc)i.com/supernovasearc)i.Fitml 

Fittp://ita)iansupemovae.org/ 

FittpV/www.weizmann.ac.il/astropFiysics/wiserep/ 
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TABLE 1 

Supernova Sample 


Name 

a(J2000) 

(deg) 

5(J2000) 

(deg) 

Z 

DM 

(mag) 

E(B-V) 

(mag) 

.3T«,pealc 

(mag) 

to 

(day) 

^peak 

(day) 

FAP 

DP 

ref. period 

PTF 09dxv 

347.144705 

+ 18.937131 

0.0322 

35.77 

0.145 

-18.0* 

55079 

55094 

0.000 


>56041 

SN2009nf, PTF09gyp 

029.736494 

-07.282473 

0.046 

36.57 

0.027 

-17.6 

55138 

55151 

0.025 


>55835 

PTF 09hnq 

345.470095 

+ 14.413534 

0.027 

35.38 

0.108 

-17.5* 

55157 

55168 

0.002 


>56449 

PTF 09ism 

176.149461 

+ 10.212143 

0.03 

35.61 

0.075 

-17.4 

55194 

55200 

0.000 


>56328 

PTF lOfqg 

190.457745 

+ 11.591142 

0.0278 

35.44 

0.032 

-16.6 

55302 

55323 

0.126 


>56712 

PTF lOqil 

347.470125 

+ 13.132566 

0.0396 

36.23 

0.076 

-17.0 

55413 

55427 

0.027 


>56503 

PTF lOtzh 

257.305622 

+41.755139 

0.034 

35.89 

0.025 

-16.3* 

55437 

55459 

0.002 

y 

>56360 

PTF lOxfl 

034.955874 

+ 15.295001 

0.05 

36.76 

0.105 

-18.2* 

55456 

55470 

0.030 


>56533 

SN201 Idh, PTF 1 leon 

202.521152 

+47.169782 

0.001683 

29.45 

0.035 

-17.1 

55713 

55732 

0.050 

y 

>56662 

SN2011hg, PTFllpdj 

347.953215 

+31.016708 

0.0236 

35.08 

0.075 


55856 


0.000 


>56735 

PTF llqju 

213.862526 

+36.408608 

0.0282 

35.48 

0.089 

-17.5 

55880 

55902 

0.007 


>56324 

SN2012P,PTF12os 

224.996177 

+01.890051 

0.004533 

32.04 

0.051 

-16.3 

55931 

55952 

0.006 


>56710 

PTF 12fxj 

027.004256 

+35.708841 

0.0148 

34.06 

0.056 

-17.4 

55089 

55102 

0.054 


<55127 

SN2012ey, PTF12iqw 

036.210492 

+ 16.181395 

0.027 

35.38 

0.179 

-15.8* 

56183 

56198 

0.000 


<55503 

PTF 12jaa 

338.364087 

+00.740223 

0.0237 

35.09 

0.087 

-15.6* 

56190 

56193 

0.000 

y 

<55896 

PTF13nu 

183.267741 

+32.613898 

0.026 

35.30 

0.013 

-17.7 

56356 

56376 

0.355 


<55896 

SN2013bb, PTF Baby 

213.058377 

+ 15.842080 

0.01755 

34.43 

0.015 

-16.5 

56369 

56394 

0.217 


>56766 

PTFBajn 

129.528692 

+66.526227 

0.03 

35.62 

0.049 

-17.3* 

56385 

56404 

0.077 

y 

<55183 












& <56325 

SN2013CU, PTF Bast 

218.495242 

+40.239672 

0.0258 

35.28 

0.012 

-18.5 

56414 

56426 

0.003 


<55042 

PTFBebs 

140.284454 

+49.592651 

0.027 

35.38 

0.022 

-16.9* 

56606 

56622 

0.003 


>56780 

SN 201 lef 

352.737583 

+ 15.490083 

0.0134 

33.84 

0.067 


55760 


0.095 


>56406 

SN2012an 

261.042625 

+59.001916 

0.0111 

33.43 

0.030 


55978 


0.000 


>56778 

SN2012CS 

232.990208 

+68.245222 

0.0218 

33.56 

0.027 


56053 


0.041 


<55052 

LSQ 12fwb 

006.355083 

+06.707194 

0.03 

35.62 

0.024 


56233 


0.114 


<55430 

LSQ 12htu 

152.904625 

-07.386555 

0.04 

36.26 

0.034 


56282 


0.006 


>56712 

SN20Bdf 

186.622208 

+31.227305 

0.002388 

31.23 

0.019 

-16.8 

56447 

56470 

0.965 

y 

<56388 

PSl-14od 

050.275958 

-07.282611 

0.02 

34.72 

0.059 


56713 


0.000 


<55089 


Note. — The SN sample. a(J2000) and 5(J2000) are the J2000.0 right ascension and declination, respectively; z is the SN redshift obtained from 
spectroscopy. DM, the distance modulus, is derived from the redshift with Hq — 69.33 kms“^ Mpc“^, Q.m — 0.24, and f2A=0.71 JHinshaw et al.l2QT^ . The only 
exceptions are the three closest SNe and SN2Q12c s , whe re redshift-independent distance measurements of the host galaxies are available on NED. EiB — V) 
is the Galactic extinction taken from ISchlegel et aT] j 19981) . jW^ peak is the absolute 7?-band magnitude of the brightest detection; asterisks indicate that the peak 
is not well observed and the SN might be considerably brighter, to is the MJD of the approximate explosion date estimated by picking a date between the last 
nondetection and the first detection; thus, for some SNe the uncertainty in to can be many days. ?peak is the MJD of the brightest detection, where only the second 
peak is considered in case of a double-peaked light curve. FAP, the fal se-a larm probability, is the probability of detecting a false precursor candidate by coadding 
images in 15-day bins as estimated using the bootstrap method (see ED- The penultimate column (DP) indicates whether two peaks are observed in the light 
curve, and the last column specifies the chosen reference period. Below the solid line we list SNe added to our sample from the literature. 

References: 

PTE 09hnq, PTF lOfqg, PTE lO qrl, PTE lOtzh, PTE lOxfi, PTE llqju, PTE 12jaa, PTE 13nu, PTE 13ajn, and PTE 13ebs: reported here for the first time. 

PTE 09dxv: lArcavi et^ 120101) . 

SN2009nf (PTF09gvp):IDrake etalJd2009bD : l^avi et al.l (2013). 

PTF 09ism: l Arcavi et al.H2O10l) ! 

SN201 Idh (PTF 1 leon): lArcavi et al.l (Ml;IGriga et al.l (Ml. _ 

SN 201 Ihg (PTF 1 lpdj)J^^^^^T&_Mazzoi0U20n^^]GaF%m^tayj2(ni^JMmon&^erliny <201 lf) : lTomasella et alj <20111) . 

SN2012P (PTF 12os): lArcavi et al.l <2012l):IBorsato & Nascimbenil <2012l) : inimai et al.l <20121) . 

SN 2013bb (PTFJ3ab^3S owertmie^ n<2013D:IElias-Rosa et alj <20131) . 

SN201 lef: IBlanchard et al.l <20TTl )i lPan'enreral I <201 ll i^ 

SN 2012an: IChen et alJ j20i3) : IJha et al.l iM^ : lNe^on & Pucketll (20T3). 

SN2012cs: fRich et alj<20m 
onetalTlB 


SN 2012ey: [Howerto^ 


<20igi:ITuratto et al.l f20Tl). 


,_ DC _. 

LSQ 12htu: ILe Guillou ef^ (M2h . 

SN2013cu: lGal-Yam et alj <20141) . _ 

SN2013df: lQaEatt^et^<201^ dVan Dvk et aT] <2013ll^(?T^ ; IMorales-Garoft'olo et al.l <2014l) : IBen-Aini et al.l <20151) . 
PSl-14od: ICMnDbeiTetalTKo 1 41) . 


by visiting selected fields of the sky regularly over the dura¬ 
tion of se veral months . The data-reduction procedure is de- 
scr ibed bviLaher et al.l (201^ and the photometric calibration 
by lOfek et al.l (l2012allb ). We only use observations in the R 
band and neglect the < 15% of data in the g band to simplify 
this analysis. For every SN a reference period containing at 
least 20 observations is chosen. If possible, we use data well 
after the SN has faded to construct the reference image, but 
when no such data is available we instead resort to the oldest 
pre-explosion images. For PTF 13ajn, a considerable number 
of observations were acquired with two CCDs, and we hence 


define two different reference periods. All reference periods 
are listed in the last column of Table [T| 

For each SN the reference image is created by coadding 
observations within the reference period. This image is sub¬ 
tracted from all science images and the (positive or nega¬ 
tive) fiux residuals at the SN location are measured using 
forced point-spread function (PSF) photometry. A correc tion 
for Galactic extinctio n according to lSchlegel et al.l (Il998l) and 
ICardelli et al.l (Il989ll is applied to all fiuxes and magnitudes 
in this paper. 

The observations and their fiux residuals are listed in 
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TABLE 2 

PTF Observations 


Name 

MJD-fo 

(day) 

MJD 

(day) 

niPTF.R 

(mag) 

mpTF,^ En' 

(mag) 

Lim Mag 
(mag) 

Flux 

(counts) 

Flux En- 
(counts) 

PTF 09dxv 

-47.876 

55031.304 

25.34 

31.38 

21.09 

2.7 

76.9 

PTF 09dxv 

-47.805 

55031.375 

22.72 

1.24 

21.39 

51.1 

58.4 

PTF 09dxv 

-45.711 

55033.469 

78.65 

- 

21.38 

-9.0 

59.2 

PTF 09dxv 

-45.705 

55033.475 

21.69 

0.49 

21.37 

133.2 

59.7 

PTF 09dxv 

-43.762 

55035.418 

79.45 

■ 

21.46 

-18.8 

54.9 


Note. — Flux residuals in the pre-explosion light curve of all SNe in our sample. Magnitudes are calculated as “asinh magnitudes” ILupton et alJI 199^ . 
and have a meaning only when smaller than the limiting magnitude. The limiting magnitude here is at the 3f7 level. This table is published in its entirety in the 
electronic version of ApJ. A portion of the full table is shown here for guidance regarding its form and content. 


TABLE 3 

Log oe spectroscopic observations 


Name 

Telescope 

Instmment 

date 

PTF 09dxv 

Keck 1 

LRIS 

2009-10-22 

SN 2009nf 

Keck 1 

LRIS 

2009-11-11 

PTF 09hnq 

HET 

LRS 

2009-11-25 

PTF09ism 

P200 

DBSP 

2010-01-09 

PTF lOfqg 

P200 

DBSP 

2010-05-07 

PTFlOqrl 

Keck 1 

LRIS 

2010-09-05 

PTF lOtzh 

Keck 2 

DEIMOS 

2010-10-12 

PTF lOxfl 

P200 

DBSP 

2010-10-17 

SN 201 Idh 

Lick 3-m 

KAST 

2011-07-06 

SN 201 Ihg 

Lick 3-m 

KAST 

2011-12-19 

PTFllqju 

Keck 1 

LRIS 

2012-02-20 

SN2012P 

HET 

LRS 

2012-01-14 

PTF12fxj 

Lick 3-m 

KAST 

2012-07-11 

SN 2012ey 

Keck 2 

DEIMOS 

2012-10-14 

PTF12jaa 

Keck 2 

DEIMOS 

2012-10-14 

PTF13nu 

P200 

DBSP 

2013-04-13 

SN 2013bb 

Keck 1 

LRIS 

2013-04-09 

PTF13ajn 

P200 

DBSP 

2013-05-02 

SN 2013CU 

P200 

DBSP 

2013-06-02 

PTF13ebs 

Keck 1 

LRIS 

2013-12-02 


Note. — The spectra are shown in Figures [SIIOI 


Table |2] They were obtained up to 3.5 years prior to the SN 
explosion. Figures[T}{3display the pre-explosion light curves 
in 15-day bins. Filled circles correspond to bins containing 
six or more observations, while open circles represent bins 
with fewer data points. Crosses around zero mark th e es- 
timated 5(7 noise level, whose calculation is described in il4.1l 


4. PRECURSOR SEARCHES 

The precursor sea rch follows to a large extent the methods 
of lOfek et al.l (120141) . In a first search we look for precursor 
events in individual images. In a second search channel ob¬ 
servations are combined once in 15 and once in 60-day bins 
to gain sensitivity to faint precursor events. An advantage of 
the binned search is that it allows to estimate the noise level 
from the data itself which is more precise than the formal er¬ 
rors calculated by the image subtraction pipeline. The only 
drawback of binning the data is that a minimal precursor du¬ 
ration has to be assumed and we loose sensitivity to shorter 
eruptions. 

Both search methods a re p resented in the 114.11 and the 
results are described in 114.21 Add ition al tests and cross 
checks on the data are explained in 114.31 This search uses 
tools ava ilable in the MATLAB astronomy and astrophysics 
package (IOfekll20r^ . 


4.1. Search Methods 

In the unbinned search we have to rely on the statistical 
errors estimated by the image subtraction pipeline. These un¬ 
certainties are based on Poisson noise and several additional 
uncertainties are neglected. Those neglected error sources in¬ 
clude errors on the fitted PSP and correlations between neigh¬ 
boring pixels induced by s moothing images before the sub¬ 
tractio n with a kernel (e.g., lAlard & LuptonI 120101 : iBramichl 
|200^. Therefore, the calculated errors often underestimate 
the fluctuations in the data and we here scale them up to ob¬ 
tain a more realistic error estimate. 

We calculate the standard deviation of the complete pre¬ 
explosion light curve using the errors from the image subtrac¬ 
tion pipeline and compare it to the act ual scatter in the data 
estimated using the bootstrap method (lEfronl[T982h . When 
applying the bootstrap method we resample the pre-explosion 
light curve 1000 times by randomnly assigning flux residuals 
to the observation times. Plux residuals can be drawn several 
times. We then calculate the average flux of every simulated 
light curve. The average fluxes are distributed according to 
a normal distribution whose standard deviation is the boot¬ 
strap error on the average flux. When using the bootstrap 
technique, we implicitly assume that the pre-explosion light 
curves are dominated by statistical fluctuations, and that even¬ 
tual precursors have a negligible influence on the light curves’ 
standard deviation. We find that the bootstrap error is always 
larger: Por many SNe the difference is approximately a factor 
of two. The errors on the individual flux residuals are scaled 
up accordingly and we require a 5(7 deviation above zero for 
the detection of a precursor. 

We caution that these flux residuals (even in the absense of 
a precursor) do not follow a Gaussian distribution. Outliers 
can be caused in a number of ways, such as by atmospheric 
conditions, cosmic rays hitting the detector, or imperfect im¬ 
age subtractions. Large deviations are rare, however, and they 
can be identified by inspecting the image or verifying that a 
precursor candidate is detected in several subsequent images. 

In the second search channel the sensitivity of the search is 
increased by coadding flux residuals in time bins. Our coad¬ 
dition method preserves more information relative to simple 
image coaddition and has com mon grounds with the op timal 
coaddit ion method descri bed bv IZackav & Ofel3 (l2015h . Pol- 
lowing lOfek et aTI ((2013), we choose a bin size of 15 days, 
which means that precursor events having shorter durations 
might be missed even if their luminosity is above the quoted 
sensitivity. In addition we repeat the analysis using 60-day 
bins to search for long-lasting faint precursors. 

Por the binned searches bootstrap errors are calculated for 
individual bins, which is more reliable than scaling up the 
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Fig. 1.— The pre-explosion light curves of the first ten SNe (names are listed in each panel; compare Table[T). Shown are flux residuals (relative to the reference 
image given in counts per 60 s exposure) coadde d in 15-day bins pr ior to the SN explosion. The right-hand axis shows apparent mpxp.R-band magnitudes. All 
measurements are in the PTF magnitude system lOfek et al]l2012al lbh. The zero point of the flux residuals is 27. Filled circles represent 15-day bins with six 
or more observations and open circles mark bins h aving less da ta. The plus signs sho w the 5 O' noise level, which is estimated based on scaled Poisson errors 
for the open circles and with the bootstrap method jEfron|[T982ll for filled circles (see 114. ll for a detailed description). If observations are found in consecutive 
bins, the plus signs are connected by a solid line. A few very large error regions are outside the plotted range to improve readability. Figures|2]and|3]show the 
pre-explosion light curves of the remaining SNe. 


Poisson errors estimated by the image subtraction pipeline. 
To get a sound error estimate with the bootstrap method, a bin 
must contain several entries. Therefore, in the following only 
bins with at least six observations are considered. 

As before, we require a 5(7 deviation for the detection of 
a precursor, with the difference that the noise level for ev¬ 
ery bin is estimated with the bootstrap method. To calculate 
a false-alarm probability we randomly combine observation 
times and flux residuals from the entire pre-explosion light 
curve 1000 times, bin them in 15-day bins, and look for pre¬ 
cursors in these scrambled light curves. The probabilities for 


detecting one or several false precursor candidates per SN are 
listed in Table [T] This test, however, only leads to a valid re¬ 
sult if the pre-explosion light curve is dominated by statistical 
fluctuations and its mean flux is consistent with zero. 

The binned pre-explosion light curves are shown in Figures 
[T]-[3] In addition, bins containing fewer than six observations 
are marked by open circles. For these bins, the standard 
deviation of the mean flux is calculated based on the scaled 
Poisson errors. 
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Fig. 2.— Figure[T] continued. 

4.2. Precursor Candidates 

When applying the precursor searches described in 114. II to 
the SNe in our sample, no precursor candidates are found in 
the unbinned search or when using 15-day bins, however, one 
precursor candidate is discovered when binning observations 
in 60-day bins. 

In the unbinned pre-explosion light curves (not shown; 
flux residuals given in Table |2]l, no single observation has a 
flux above the 5(7 noise level estimated based on the scaled 
Poisson errors. From the nondetection of events above this 
threshold, we can infer an upper limit on the false-alarm rate 
for this sea rch channel. U sing the Poisson single-sided up¬ 
per limits (lGehrelslll986h . the expected number of precur¬ 
sors is smaller than 2.3 at the 90% confidence level when 
no event is observed. Since we evaluate in total 3152 pre¬ 


explosion observations, the false-alarm probability is smaller 
than 2.3/3152 k, 0.0007 per observation. 

We find that even when lowering the threshold of this 
search to 3(7, only one single flux residual reaches the 
threshold. The most significant subthreshold event is a 
4.2(7 deviation in the light curve of SN2012P, 654 days 
before the SN explosion. A second observation with similar 
limiting magnitude only two hours earlier, however, does 
not yield a detection. The two following observations 3 
days later are just as deep, but neither show signs of a 
flare. Given the number of trials, this event is most likely 
spurious. We nevertheless conservatively maintain the 5(7 
threshold when calculating limits on the precursor rate in @ 
Potential precursor events below the quoted flux threshold 
of our search do not affect the validity of the calculated limits. 
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Fig. 3.— Figures[T]and|2] continued. Flere we show the PTF pre-explosion light curves of SNe lib found in the literature. 


We also do not find precursor candidates above the 5(7 
noise level when using 15-day bins. The only exception is the 
pre-explosion light curve of SN2011dh, the nearest object 
in our sample. Its progenitor star is detected in coadded 
images, and in 17 of the bins the noise level is low enough 
to allow a signihcant detection. The average flux of the 
pre-explosion light curve of SN2011dh is 17(7 abo ve z ero 
and the progenitor is hence clearly detected. In il5.1l we 
account for the progenitor’s average flux and all precursor 
candidates vanish, as can be seen in Figurej^ 

When increasing the bin size to 60 days one precursor can¬ 
didate is detected marginally above the 5(7 threshold in the 
pre-exposion light curve of SN2012cs. Figure |4| shows the 
binned light curve as well as the unbinned observations. For 
60-day bins the noise level is low enough to show a flux ex¬ 
cess approximately 680 days before the SN explosion. The 16 
observations in this bin are nearly all at the same flux level. 
The excess is, hence, not caused by individual observations 
but many points contribute to it as it is expected from a long 
lasting precursor event. The median flux in this bin is 62 
counts which corresponds to a magnitude of mpTp^s = 22.51. 
SN 2012cs has a redshift of 0.0218. The absolute magnitude 


of the possible precursor is Mpxp = —11.0 and its luminos¬ 
ity is Lprec. = 7 X 10^® ergs^* (no bolometric correction ap¬ 
plied). A duration of about 60 days seems consistent with 
the unbinned light curve and the radiated energy amounts to 
4 X 10^® erg during this period. If real, the event is hence 
slightly fainter th an precursors detected so far (compare e.g., 
lOfek et aTll2014ft and its radiated energy is similar to the least 
energetic detected precursors. 

Since the event is barely above the signihcance threshold 
of this search it might be a false detection. The false-alarm 
probability estimated for SN 2012cs and a bin size of 60 days 
is 1.8%. Its calculation is, however, based on the assump¬ 
tion that the pre-explosion light curve is dominated by statis¬ 
tical fluctuations. If the event is real, a considerable fraction 
of the observations (16 out of 52) are systematically shifted 
upwards. A real precursor can thus enhance the false-alarm 
probability. We conclude that the detected candidate is con¬ 
sistent with an astrophysical precursor event. However, due to 
the marginal detection we cannot decide whether or not it is 
real. No further precursor candidates are found. We focus in 
the following on the nondetection of precursors using 15-day 
bins, which allows us to constrain the rate of bright precur¬ 
sors. The detected possible precursor is below the sensitivity 
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Fig. 4.— A section of the pre-explosion light curve of SN2012cs. The 
gray error bars show the scaled up 1(7 Poisson errors on the individual ob¬ 
servations. The blue lines indicate the bin position as well as the 5(7 noise 
level for the binned light curves. Open dots represent bins containing fewer 
than 6 obseiwations for which the noise level is estimated based on the scaled 
Poisson en'ors. Filled dots correspond to bins with at least 6 observations 
and the bootstrap method was used to estimate the noise level. A detection, 
marginally above the 5(7 threshold, is obtained around day —680. The false 
alarm probability of SN 2012cs is 1.8% when using 60 day bins. 

threshold of this analysis and hence does not affect the limits. 
The pre-etmlosion light curves of all 27 SNe Ilb are shown in 
Figures [T]-[3 and the upper edge of the error band corresponds 
to the sensitivity of the search with a bin size of 15 days. 

4.3. Additional Tests 

In addition to the tests described in il4.1l we check whether 
the pre-explosion light curves are consistent with an average 
flux of zero, and we look for autocorrelations in the unbinned 
light curves. Both tests, described below, do not reveal hints 
for detections or raise concerns about the data quality. 

Initially, the complete pre-explosion light curves of several 
SNe were inconsistent with zero when calculating its aver¬ 
age flux and its bootstrap error. In some cases the SN was 
still present in the chosen reference images, and we instead 
resorted to a pre-explosion reference. For PTF lOqrl, we use 
pre-explosion images since the SN location is at the very edge 
of the CCD in the post-explosion images, affecting the image 
subtraction. Also, the average pre-explosion fluxes of several 
sources show deviations marginally below the 3(7 threshold. 
This may affect the false-alarm probabilities of the binned 
search (shown in Table [T]i. We verify, however, that no ad¬ 
ditional precursor candidates are found when subtracting the 
mean pre-explosion flux from all of the residuals. 

Furthermore, we look for autocorrelations in the unbinned 
pre-explosion light curves where we calculate the correla¬ 
tion between an observation and the five following observa¬ 
tions. The only deviation above the 3(7 threshold is found for 
PTF lOqrl with a lag of two observations. This might, how¬ 
ever, be caused by chance. The autocorrelation is calculated 
for all 27 SNe with five different lags each, which means that 
at the 3(7 level 0.36 false detections are expected. The proba¬ 
bility to detect one or several events above the 3(7 threshold is 
therefore 30%. Moreover, when looking for autocorrelations 
in time bins of three days instead of indexing the observations, 
the correlation vanishes. PTF lOqrl is located at a redshift of 
0.05 and is hence is too far away for the detection of the pro¬ 
genitor star. 


We conclude that no indications for precursor detections are 
found in these additional tests. 

5. PROGENITOR DETECTIONS AND VARIABILITY 

Eor the three nearest SNe in our sample, here we explore 
the possibility of detecting the progenitor star and calculating 
limits on its variability. 


5.1. SNlOlldh 

SN 201 Idh was discovered on 2011 May 31 in the nearby 
galaxy M51 at a distance of 7.8 Mpc, less than 15 hours 
after a nondet ection down to a limiting magnitude of 
mpjp g = 21.4 (lArcavi et alJ 1201 lb . The progenitor star is 
clearly detected both by the Hubble Space Telescope ( HST ; 
iLi et all 1201 U iMaund et all 1201 U iVan Dyk et all 1201 ih and 
in gro und-based obseryation s (e.g., the Large Binocular Tele- 
scoDe.l^z czygiel et alJl2012l : the Northern Optical Telescope, 
lErgon et^ 201^ ; and PTF). The progenitor is a yellow su- 
pergiant (IVan Dyk et al.ll2013l) . Moreoyer, a flux excess has 
been measured in the fading SN light curye in blue bands 
and has been attributed to the presence of a binary compan¬ 
ion star (iFolatelli et al.ll2()T4li . The companion presumably 
hasMF 225 w = —5.11 ±0.29 mag andMpsaew = —4.66 ±0.29 
mag, consistent with a B star haying IOMq < M < I 6 M 0 
(IFolatelli et akllMl) . 

There are 373 PTF R-band obseryations that span the last 
two years before the explosion of SN 201 Idh. Our reference 
is based on obseryations obtained starting 949 days after the 
explosion, when the SN has faded away. The progenitor star 
is hence not present in the reference image, so its flux is not 
subtracted from the pre-explosion light curye. All flux resid¬ 
uals are listed in Table |2l 

The progenitor is not detected in indiyidual images; how- 
eyer, in the binned light curye shown in Figure 2] the SN 
is significantly detected in many bins. The ayerage un¬ 
absorbed flux of the progenitor is 139 ± 10 counts, cor¬ 
responding to mpjp]( = 21.82 ±0.08 mag, where the un¬ 
certainty is the 1(7 bootstrap error. During the HST ob¬ 
seryations in 2005, the progenitor star was detected in the 
F658N band with a magnitude of mpessN = 21.39 ±0.02 
(IVan Dyk et al.l 1201 1[) . corresponding to a flux of (4.94 ± 

0.09) X 10^** erg s^^ cm^^ A ^ at an effectiye wayelength of 
6579.5 A (flux taken from iBenyenuto et al.ll2013h . This yalue 
is comparable to our measurement, translating into a flux of 

(5.54 ±0.35) X 10^'^ erg s^^ cm^^ A * at the same effectiye 
wayelength. 

To search for precursors on top of the progenitor flux, the 
5(7 error band is centered around the star’s ayerage flux. No 
single or binned obseryations exceed this error band, and 
hence no precursors are found. One bin around day 215 before 
the explosion has an ayerage flux of more than 600 counts. 
It howeyer only contains two flux residuals which both haye 
large errors. We hence conclude that this point is caused by 
two low-quality obseryations. Its deyiation from zero is less 
than 2 ( 7 . 

Moreoyer, we search for periodic yariability by calculating 
the power spectrum of the pre-explosion flux residuals. No 
significant yariability was found at any time scale either in 
th e binned or unb i nned d ata. 

ISzczygiel et al.l (1201 2 l) analyze pre-explosion obseryations 
taken with the Large Binocular Telescope within three years 
before the SN explosion. With data at flye different epochs 
they suggest a continuous fading at a rate of 0.039 ± 
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Fig. 5.— The binned pre-explosion light curve of SN2011dh. The ref¬ 
erence images are chosen to be after the explosion, hence the flux of the 
progenitor has not been subtracted. No observations deviate by more than 5(7 
from the mean flux. 

0.006 mag yr^* in the R band. However, in a sample of stars 
of comparable luminosity, as large or larger dimmings are ob¬ 
served for 5% of the stars, and about 16% of the stars have 
a photometric root-mean -square value larger th an that of the 
progenitor of SN 201 Idh (ISzczvgiel et alJ2012h . Therefore, it 
is possible that this variability is not real. A fading in the last 
century before the SN explosion is expected as the envelope 
slowly responds to the rapid changes in the stellar core during 
the last stages of nuclear burning. The decline rates predicted 
by simulations are, however , typically 100 times s maller than 
the measured dimming rate (ISzczvgiel et alJl^l2^ . 

The measured dimming of 0.039 mag yr^* would corre¬ 
spond to a change of only 9 counts between our hrst and last 
observation, comparable to the uncertainty in the average 
flux. To check whether our observations do favor a dimming, 
the binned pre-explosion light curve is htted both with a 
constant flux and with a dimming at the measured rate. 
The results are very similar, with X^/dof = 35.4/39 for 
the constant model compared to x}jdof = 34.4/39 for the 
fading light curve. With g a 1, we cannot conhrm or 
eliminate the fading hypothesis of ISzczvgiel et al.l (120121) . 

5.2. SN2013df 

The second-nearest object in our sample is SN2013df at 
a distance of 18.2 Mpc. No attempt can be made to detect 
the progenitor since we do not have potential reference 
images after the SN has faded. The progenitor, a^OOR© 
cool supergiant, was identihed by IVan Dvk et al.l (l2014h in 
pre-explosion HST images with an absolute magnitude of 
Mhst.v ~ — 6 . 8 . 


5.3. SN2012P 

At the position of SN 2012P, a potential progenitor star was 
detected in archival H ST observations from 2005 March 10 
(I Van Dvk et al.ll2()I^ . The source has an absolute magnitude 
of Mhst.v ~ —9.1, and its luminosity and colors are consistent 
with a very luminous supergiant star or a star cluster. 

Interpolating the V and /-band magnitudes reported by 
I Van Dvk et al.) (120121) . a PTF /?-band magnitude of mpTPR = 


22.74 is expected, corresponding to a signal of 50 counts 
per observation. The weighted average flux of the 308 pre¬ 
explosion observations is, however, —17 ± 10 counts per ob¬ 
servation (1 (7 boo tstrap error). Hence, we conclude that the 
source reported by I Van Dvk et al.l (120 12^ is still present in the 
reference image, which consists of observations more than 
779 days after the SN explosion. The source close to the SN 
po sition could be either a compact star cluster as mentioned 
bv IVan Dvk et al.l (|20T3), a luminous binary companion, or an 
unrelated object. 

Based on the PTF pre-explosion observations in the last 
three years before the explosion of SN 2012P, we can exclude 
precursor explosions with an absolute /?-band magnitude of 
44pTF,fi <—11 (T>8x lO^^ergs^^ without a bolometric 
correction) during 30% of the time. 

6 . CONTROL TIME AND PRECURSOR RATES 

Eollowing lOfek et'H] (120141) . we combine the results of the 
individual SNe and calculate an upper limit on the precursor 
rate of SNe lib. Hereby, we implicitly assume that the pro¬ 
genitors of SNe lib form a uniform population and that our 
SN sample provides a good representation of this population. 
Both presumptions are likely not fulhlled and the resulting 
rate estimate should hence be regarded with caution. 

To address the inhomogeneity of the progenitor population, 
we additionally calculate an upper limit only for SNe lib with 
observed double-peaked light curv es since this is a stro ng in¬ 
dication for an extended envelope (iNakar & Piroll201^ . 

Moreover, faint SNe lib are likely underrepresented in this 
sample; however, our sample does cover a wide range of peak 
luminosities and light-curve shapes. This analysis is thus not 
restricted to any certain subgroup of SNe lib. The advantage 
of a rate estimate is that we potentially gain information on 
the SNe lib as a class and we are able to quantify the results 
of this analysis. 

To calculate a precursor rate for all SNe lib, we estimate 
the control time — the time during which a precursor above 
a certain absolute magnitude can be detected. The following 
calculations are described for a bin size of 15 days. When us¬ 
ing 60-day bins all calculations and results are similar except 
for the detection of the faint precursor candidate. All resulting 
upper limits or precursor rates are summarized in Table |5] for 
both bin sizes. 

We adopt a minimal precursor duration of 15 days and bin 
the observations accordingly. Eor bins with six or more obser¬ 
vations, the 1(7 noise level is calculated using the bootstrap 
method and is then multiplied by a factor of 5 to estimate 
a 5(7 level. Eor bins with fewer observations, we determine 
the standard deviation based on the scaled Poi sson errors esti¬ 
mated by the image-subtraction pipeline (see il4.11 l. The upper 
edge of the noise level corresponds to the 5(7 limiting magni¬ 
tude per bin. The control time is formally given as 

Niiii Asn 

Learch(44i?) = h,search 

if Mi g < Mg 
otherwise, 

where Mg is the absolute magnitude at which the control time 
is calculated and Mi g is the 5(7 limiting magnitude in bin i. 
The sum runs over all Abin light-curve bins per SN and all Nsn 
SNe in this sample. The mean observation time of the bin 


;=1 


with 6,search ) — 


15 days, 
0 days, 
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TABLE 4 

Precursor Search Control Time 


Name 

At 

(day) 

»IPTF,K 

(mag) 

Abs. A/pxF,y? 
(mag) 

^meas 

PTF 09dxv 

-40.0 

21.49 

-14.28 

20 

PTF 09dxv 

-25.5 

20.90 

-14.87 

23 

PTF 09dxv 

-14.4 

22.38 

-13.39 

14 

SN2009nf 

-63.4 

20.61 

-15.96 

4 

SN2009nf 

-Al.l 

22.87 

-13.70 

8 

SN2009nf 

-28.6 

20.08 

-16.49 

3 

SN2009nf 

-18.0 

20.41 

-16.16 

3 

PTF 09hnq 

-117.5 

20.50 

-14.88 

9 


Note. — Listed are the SN name, the mean observation time of the 15- 
day bin relative to the explosion date, the apparent and absolute magnitude 
down to which precursors can be excluded, and the number of observations 
per bin. The limiting magnitudes are at the 5 O' level estimated from Poisson 
errors for b ins w ith less than 6 observations and with the bootstrap method 
otherwise ( i]4.U . They correspond to the upper + signs in FiguresHUl] and 
their cumulative distribution is shown in Figure The full table is available 
in the online version. 

relative to the explosion date, the limiting magnitude, and the 
number of observations per bin are listed in Table|4](unbinned 
observations given in Table|2]i. _ 

From the search results described in 114.21 and the light 
curves sho'wn in Figures[T]-[3] we know that during the control 
time no precursors were detected above the estimated noise 
level when using 15-day bins. The control time fseaich for the 
complete sample as a function of the precursor magnitude is 
depicted in Figure |6] where fsearch is defined as the bin du¬ 
ration of 15 days multiplied by the number of bins in which 
a precursor of the respective magnitude would have been de¬ 
tected. For example, at an absolute magnitude of —14 and 
considering the full time range of 3.5 years prior to the SN 
explosion (thick black line in Figure |6]), a control time of 9.2 
years has been accumulated; in our sample there are 223 15- 
day bins in which precursors having an absolute magnitude of 
— 14 or less can be excluded. 

According to the amount of data and the SN distance, the 
objects in our sample contribute differently to the magnitude- 
dependent control time. The five SNe with the longest con¬ 
trol time (PTFllqju, SN2011dh, PTF12fxj, SN2012P, and 
SN2011ef) account for close to 50% of the total time cov¬ 
ered by observations. Especially at large absolute magnitudes, 
the control-time distribution is dominated by the two nearby 
events, SN2011dh and SN2012P. In Figure |6] we show how 
the distribution changes when removing SN2011dh, which 
owing to its small distance has a large contribution to the con¬ 
trol time in the absolute magnitude range —7 to —10. In a 
similar way, the second bump at magnitudes from —9 to —12 
is due to the observations of SN 2012P. At lower magnitudes 
many other SNe start to contribute to the distribution, so the 
results are less dominated by individual objects. 


For searches in which no precursors are detected, we can set 
a 90 % upper limit o n the expectation number of precursors of 
2.3 (lGehrelslll986t one-sided Poisson statistics). The limits 
on the average precursor rate per SN per year, Rprec, and on 
the average number of precursors per SN, Aprec, are given by 

f?prec<^ and Ap,ec<2.3^^5^, (2) 

^search ^search 

respectively, where fwindow is the duration of the considered 



Abs. Mp.^p p, [mag] 


Fig. 6. — The absolute-magnitude dependent control time for the complete 
sample, during which the SN locations were observed and no precursors were 
detected above the given magnitude. The various curves are for precursors 
taking place at different time ranges prior to the SN explosion (see legend). 
The curve normalization and shape are determined by the amount of data and 
the SN distances. SN 201 Idh accounts for most observations at the brightest 
magnitudes and the dotted line shows the effect of removing it from the sam¬ 
ple. Additional curves display the control time for shorter periods before the 
explosion. The limiting magnitudes per bin for all SNe are listed in TableO 

time interval before the explosion. The limit on Aprec is hence 
purely determined by the ratio fwindow/Search . the amount of 
time during which precursors can be excluded compared to 
the time window during which the search is conducted. We 
note that Rprec is an average quantity since the bins of fsearch 
are not necessarily distributed homogenously within fwindow- 

The resul ting upper limit on the precursor rate is shown in 
Figure [7(a)l and Table|5]compares the upper limit on the num¬ 
ber of precursors pe r SN to the numb er of precursors mea¬ 
sured for SNe Iln bv lOfek et al.l (l2014l) . 

The upper limit on the number of precursors for the longest 
time interval of 3.5 years prior to the SN explosion is smaller 
than unity. This means that at the 90% confidence level, 
not all SNe Ilb can exhibit precursors as bright as —14 mag, 
which is a typical /?-band luminosity of precursor s detected so 
far (e.g.. IPastorello et alJl2 007l : lOfek et alJl20T4l) . Compared 
to SNe Iln (lOfek et al.ll2014l) . the average number of precur¬ 
sors prior to SNe lib at m^nitude —14 is lower by at least 
a factor of five (see Table |5jl. In case precursors are equally 
common before SNe llb and SNe lln, they have to be at least 
as dim as an absolute magnitude of —10 for SNe llb, about 
a factor of 40 dimmer compared to the precursors observed 
prior to SNe Iln (see Tabled. 

As explained in the introduction, a double-peaked light 
curve presumably requires the presence of an extende d 
hydrogen envelope (iBersten et al.ll201^lNakar & Piroll2014!) . 
The duration of the first peak is usually only hours to few 
days and most SNe are not discovered as early. For most 
SNe in our sample we thus do not know whether or not an 
early peak is present. Ho wever, not all SNe llb have an early 
peak (lArcavi et al.ll2015]) . and there a re indications that some 
SN ll b progenitors might be compact (IChevalier & SoderbergI 
I 2 OIOI) . To account for this diversity in the progenitor popu¬ 
lation, we additionally calculate limits considering only the 
data of the five SNe llb for which two peaks are observed 
by PTF or reported in the literature. As stated in Table [T] 
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Abs. Mp^p pj [mag] 

(a) 90% upper limit on precursor rate derived from the complete SN lib sample 



Abs. Mp.pp p [mag] 

(b) 90% upper limit on precursor rate derived from five double-peaked SNe lib 


Fig. 7.— The upper limits on the average precursor rate per SN calculated based on the control time shown in Figure]^ At the lowest (brightest) magnitudes 
the limit is dominated by observations of SN 201 Idh. Considering only double-peaked SNe, the limits are weaker owing to the smaller amount of data. 


TABLES 

A^piec, AVERAGE NUMBER OE PRECURSORS PER SN 


Time interval 

SNe lib, bright precursors (—14 mag) 

SNe lib, faint precursors (—10 mag) 

SNe tin (-14 mag) 

(yr) 

15-day bins 

60-day bins 

15-day bins 

60-day bins 

15-day bins 

-3.5 to 0 

<0.86 

<0.51 

< 5.03 

, c9+3.50,+7.12 
^•^^-1.26,-1.49 

- 

-2.5 toO 

- 

- 

- 

- 

1 q-F3.8,+7.6 

4-o_2.3-3.5 

-1 toO 

<0.56 

<0.36 

<2.95 

<2.00 

^ o-F2.2,+4.4 
^•°-1.3,-2.1 

-1/3 to 0 

<0.52 

<0.29 

< 2.08 

<2.33 

ry ff-F2.0,+4.0 
"^•^-1.2,-1.8 

-2 to -1 

<0.70 

<0.47 

< 3.12 

, p,-F2.33,+4.75 
^•^^-0.84,-0.99 



Note. — Upper limit on allowed number of precursors N^rec within the respective time interval (rates from Figure [/(^ multiplied by time interval) compared 
to measured number precursors for SNe Iln <Ofek et al.ll2014h . The limits apply to precursors with a minimal duration of 15 and 60 days respectively. The 
constraints on long precursors are stronger since more observations a re coadded wh ich leads to a better sensitivity. All limits are at the 90% confidence level and 
the quoted uncertainties correspond to the 1(7 and 2(7 Poisson errors IGehrelsIl 19861) . 


this includes PTFlOtzh, SN 201 Idh, PTF12i aa, PTF13ajn, 
and SN2013df (see also lArcavi et al.l l2015l) . These SNe 
contribute only 15% of the control time of the complete 
sample, and th e limi ts are hence drastically weakened, as 
shown in Figure [7(b)| 


7. DISCUSSION 

In this section we estimate whether the explosive ejection 
of a low-mass stellar envelope in a precursor event is likely to 
be bright enought to be detectable in this analysis. We expect 
a radiatively efficient precursor eruption to radiate an amount 
of energy comparable to (or larger than) the binding energy of 
the ejected envelope, which is given by 


^bind 


^ ^env ^core 
f^core 


(3) 


where G is the gravitational constant, Menv and Mcore are re¬ 
spectively the masses of the envelope and the core, and Rcove 
is the core radius above which the envelope is located prior to 
its ejection. 

According to Nakar & Piro (2014; see also lPiro 12015h . both 
Menv and Rcoie can be estimated from the shape of the first 


peak in the optical light curve. Menv is determined by the 
time at which the bolometric light curve reaches the first peak 
(fp), and Rcore can be derived from the minimal luminosity 
Train between the two peaks. Lmin can be obtained relatively 
precisely whenever a first peak is observed. The rise time fp 
is usually less well constrained, since it can be as short as a 
day or several hours, and in most cases the data are is not 
sufficient. 

In our sample, the only double-peaked SN for which a good 
upper limit on fp is available is SN 201 Idh, where the maxi¬ 
mum of th e first peak was rea ched at most 15 hours after the 
explosion (lArcavi et alJl20111) . This limit refers to the peak 
in visible light. However, the bolometric peak probably takes 
place even earlier, henc e this limit ho l ds. Fo r SN 201 Idh, nu¬ 
merical simulations by iBersten et al.l (120121) suggest Mcore ~ 
4M0, Menv ~ 0.003 Mq, and Rmn-. » 5 x 10| * cm, c ompara- 
ble to the values estimated by iNakar & Pirol (|20T3). Using 
Eq. [2 the binding energy is Tbind ~ 8 x 10^^ erg. This re¬ 
sult is an order-of-magnitude estimate since most quantities 
entering the calculation have large uncertainties of up to a 
factor of a few. The kinetic energy required to unbind the 
envelope relates to the bolometric energy released in the pre¬ 
cursor event, Tprec, via an unknown radiation efficiency factor 
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^rad,prec — ^ ^bind (sCc lOfek 6t al.| [2 ()15l) . 

In the case of SN2011dh, within 2 years prior to the SN 
explosion, 50% of the time is covered by observations with a 
limiting magnitude of —9 or less and gaps no longer than two 
weeks. This limits the precursor luminosity during this time 
to Lprec < lO^^ergs^^ (no bolometric correction applied). If 
the envelope of this progenitor star is unbound during our ob¬ 
servations (50% probability), we can set the following limit 
on the radiative efficiency e, as a function of the precursor 
duration 5t, where the emitted luminosity is assumed to be 
constant over the precursor duration: 

5t y 

^bind ^bind \ daySy 

The unbinding event (if it ever existed) might be too faint to 
be detected if the efficiency e is low or if the precursor lasts 
for several months or longer. Another possibility is that the 
envelope is ejected prior to the time interval probed in this 
analysis. 

The limit on the precursor luminosity Lprec is close to 
the Eddington luminosity fo r a 4M(^ pro genitor. Thus, a 
continuum-driven wind (e.g. JShavivll2001h . which could in 
principle explain giant eruptions of luminous blue variable 
stars or SN Iln precursors, is likely ruled out as an expla¬ 
nation for the Type Ilb SN envelope ejection. Stellar winds 
can, however, be driven by other mechanisms such as line ab¬ 
sorption or abso rption by dust which do not require as high 
luminosities; see lLangCT(l2012^ for a review. 

If the precursor candidate detected prior to SN2012cs is 
real it released an energy of ~ 4 X 10^^ erg (see ^4■.2i . We do 
not have much information about the progenitor star. How¬ 
ever, if it is similar to the progenitor of SN 201 Idh the radi¬ 
ated energy is similar to the required energy for the unbinding 
event. It is hence possible that a stellar envelope was ejected 
in this event. 

8. SUMMARY 

We examine the pre-explosion light curves of 27 nearby 
SNe lib, searching for outbursts. One precursor candidate is 
marginally detected in a single 60-day bin in the pre-explosion 
light curve of the nearby SN 2012cs. The probability to mea¬ 
sure such an event caused by noise is 1.8% for this pre¬ 
explosion light curve. The possible precursor happened 680 
days before the SN explosion and if real its absolute R-band 
magnitude is Mji pjp = —11.0, which corresponds to a lumi¬ 
nosity of 7 X 10^® ergs^* (no bolometric correction applied) 
and a radiative energy release of 4 x 10^^ erg within the ap¬ 
proximate duration of 60 days. 

When binning the observations in 15-day bins no precursor 
eruptions are found, and we calculate a magnitude-dependent 
limit on the average precursor rate among SNe lib. Precursors 
as luminous as —14 magnitude occur on average <0.86 times 
within the final 3.5 years before the SN explosion, while in the 
last year the average number of precursors is limited to < 0.56 
at 90% confidence level. These limits are obtained under the 
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assumption that precursor eruptions last for at least 15 days. 
We conclude that bright precursor explosions, if they occur at 
all, are rare and do not happen before the explosion of every 
SN of Type lib. 

By contrast, precursors are frequent for SNe Iln (lOfek et al.l 
I20I4I) : at the 99% confidence level, 98% of all SNe Iln have 
precursor eruptions brighter than absolute magnitude —14 
within the 2.5 years befo re the explo s ion. T he precursor rate 
of SNe Iln measured by lOfek et akl (l20I4h exceeds the up¬ 
per limit for SNe lib by about a factor of five at an absolute 
magnitude of —14 or less. In addition, the precursor rate for 
SNe lib can be constrained at higher magnitudes of up to —7. 

For the nearby SN lib SN 2012P, our observations show that 
the source de tected in pre-explosio n HST observations at the 
SN position (I Van Dvk et al.ll20T2h is still present more than 
two years after the SN explosion and thus cannot be the pro¬ 
genitor star. 

The progenitor of SN 201 Idh, the closest SN in our sam¬ 
ple, is clearly detected in coadded observations, and its mean 
R-band magnitude is = 21.82 ±0.08 mag, consistent 

with archiv al HST observat i ons. T he possible slow fading 
reported by ISzczvgiel et al.l (1201 2h is below the sensitivity 
threshold of our observations. However, with 373 observa¬ 
tions over the last two years prior to the SN explosion the pro¬ 
genitor is monitored nearly constantly, and no variability or 
precursors are detected. We argue that for this progenitor star, 
the ejection of the stellar envelope in a precursor event might 
be observable except if the process is radiatively inefficient or 
lasts over several months. 
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Rest Wavelength of He features (A) 


Fig. 8.— Representative spectra of the PTF-detected SNe in our sample. All spectra are redshifted by the indicated velocity relative to the rest frame such that 
their characteristic helium absorption features are at their rest wavelengths at 5876 A, 6678 A, and 7065 A (indicated by broken lines). In addition, the black 
triangles mark the prominent Ha absorption lines. The Ha emission line is found to the right of these features and often has a flat-top profile owing to the helium 
absorption at 6678 A. Some SNe are compared to known SNe lib identifie d with SN IP ^Blondil^& T onrvll2007n and the time after peak is indi cated next to all 
spectra. All spectra are available electronically via the WISeREP webpage lYaron & Gal-Yaml2012l) . Further spectra ai‘e shown in Figures!^and llQI 
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Rest Wavelength of He features (A) 


Fig. 9.— Figure!^ continued. 
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Rest Wavelength of He features (A) 


Fig. 10. — Figure[3and[^ continued. 
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